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Abstract 

The cross section for anti-deuteron photoproduction is measured at HERA at a mean centre- 
of-mass energy of W^-^p = 200GeV in the range 0.2 < px/M < 0.7 and \y\ < 0.4, where 
M, pt and y are the mass, transverse momentum and rapidity in the laboratory frame of the 
anti-deuteron, respectively. The numbers of anti-deuterons per event are found to be sim- 
ilar in photoproduction to those in central proton-proton collisions at the CERN ISR but 
much lower than those in central Au-Au collisions at RHIC. The coalescence parameter 
B2, which characterizes the likelihood of anti-deuteron production, is measured in photo- 
production to be 0.010 ± 0.002 it 0.001, which is much higher than in Au-Au collisions 
at a similar nucleon-nucleon centre-of-mass energy. No significant production of particles 
heavier than deuterons is observed and upper limits are set on the photoproduction cross 
sections for such particles. 
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1 Introduction 



This paper describes a measurement of the rate of production of anti-deuterons in photon-proton 
collisions at HERA performed by the HI Collaboration. The measurement is of particular inter- 
est in the context of recent studies of heavy ion collisions [1]. In these collisions, the deuteron 
and anti-deuteron production rate is thought to depend on the dimensions of the collision "fire- 
ball" at the stage at which the hadrons decouple [2], i.e. when final state interactions become 
unimportant. This is the so-called "thermal freeze-out" region. 

Further, a search is performed for the photoproduction of unknown charged stable heavy 
particles in the highest energy electron-proton collisions currently accessible in the laboratory. 
The discovery of such particles would be an indication for physics beyond the Standard Model. 

The production of nuclei in particle collisions can be described in terms of the coalescence 
model. In this model [3], the cross section, a a, for the formation of an object with A nucleons 
with total energy Ea and momentum P, is related to that for the production of free nucleons in 
the same reaction, a^, with energy E^^ and momentum p = P/A, by 



where Ba is the coalescence parameter, which is inversely proportional to the source volume in 
heavy ion collisions [2], and a is the total interaction cross section of the colliding particles. 

For the measurements described here, the particles are identified through a combination of 
their specific ionisation energy loss, dE/dx, and their momenta. The numbers of anti-deuterons 
are measured in the laboratory frame in the rapidity region \y\ < 0.4 which corresponds to a 
centre-of-mass rapidity of between 1.6 and 2.4 units ^ In this range, the multiplicity distri- 
butions are on the central plateau [4] and so comparisons can reasonably be made with the 
measured numbers of anti-deuterons per event in central proton-proton collisions [5, 6]. Both 
are contrasted with data from heavy ion collisions [7, 8]. These comparisons are restricted to 
central collisions at a centre-of-mass energy greater than 50 GeV, i.e. well above the thresh- 
old for anti-deuteron production. There have been several other measurements of anti-deuteron 
production in proton-proton (pp) [9], proton-nucleus (j)A) [10], nucleus-nucleus (A A) [8] and 
electron-positron [11] collisions. These measurements are either for non-central production or 
are at a centre-of-mass energy below 50 GeV. 

2 Experimental Procedure 
2.1 The HI Detector 

Collisions of 27.6 GeV positrons with 820 GeV protons at HERA are detected in the HI detec- 
tor, which is described in detail elsewhere [12]. The components of the detector important in 

'The rapidity and pseudorapidity are defined by y = 0.5 \n[{E + Pz)/{E — p^)] and rj — ~ ln(tan ^^/2) for a 
particle with total energy E, z component of momentum p^, and polar angle 9. The -|-2;-axis (forward direction) is 
taken to be along the proton beam direction. 




(1) 
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this analysis are the small angle positron tagger, the central tracker, the backward Spaghetti-type 
calorimeter (SpaCal) and the liquid argon (LAr) calorimeter. 

The positron tagger, located at 33 m from the interaction point in the outgoing positron 
beam direction, is used to trigger on photoproduction events and to measure the energy of the 
scattered positron, from which the total photon-proton centre-of-mass energy, W^p, is deduced. 
The central track detector, surrounding the 9 cm diameter aluminium beam pipe of thickness 
1.7 mm, consists of concentric central jet drift chambers (CJCs) with inner (CJCl) and outer 
(CJC2) chambers and two additional drift chambers which measure the z coordinates of tracks. 
The pseudorapidity range covered by the central track detector is \r]\ < 1.5. The CJC has 56 
sensitive wire layers: 24 in CJCl covering radii from 20.3 cm to 45.1 cm and 32 in CJC2 
covering radii from 53.0 cm to 84.4 cm. The detector is placed inside a uniform magnetic 
field of 1.15 T, allowing measurements of the track transverse momentum to be made with a 
resolution of o-p^/pr ~ 0.009 • prfGeV] © 0.015. The specific energy loss, dE/dx, of the 
charged particles is also measured in this detector with a resolution a (dE/dx) /{dE/dx) of 
7.5% for 56 hits on a minimum ionising track. The LAr calorimeter covers the angular range 
4:° < 9 < 154° with the forward region defined to be 4° < ^ < 25°. The SpaCal calorimeter 
covers the backward region, i.e. the angular range 153° < 9 < 177.8°. The luminosity is 
measured via the well understood Bethe-Heitler process, ep — > ep7, using a photon detector at 
0° to the positron beam direction. 

2.2 Trigger Conditions, Event and Track Selection 

The measurements presented here are based on HI data taken with minimally biased triggers in 
1996 and correspond to an integrated luminosity of 5.53 ± 0.11 pb~^. Photoproduction events 
are triggered by requiring the presence of tracks in the CJC and of a scattered positron in the 
positron tagger, which ensures that the photon virtuality Q'^ < 10'^ GeV^. The following 
selection criteria are applied in order to reduce the background contamination and to ensure 
good reconstruction of the event kinematics. The selected events are required to lie within the 
interval 165 < W^p < 252 GeV (average {W^p) = 200 GeV). In this range the total acceptance 
of the positron tagger for photoproduction events is 0.46 ± 0.02. In addition, five or more tracks 
are required to be reconstructed in the CJC. These tracks must point to a common vertex with z 
coordinate within ±30 cm of the nominal interaction point. 

Candidate tracks for particle identification are selected in the range \r]\ < 1 so that they 
are well contained within the CJC. Here the track reconstruction and particle identification 
efficiencies are high. Two track selection schemes are used: the "hard" and "soft" selections. 
The minimum ionising particle (MIP) background is larger in the soft selection scheme than in 
the hard selection scheme. The hard selection is used to search for rare heavy particles, when 
it is necessary to minimise this background and to have optimum dE/dx resolution. The soft 
selection is used only for copiously produced particles such as protons and anti-protons. 

In the soft selection, the tracks are required to have at least 10 hits and to have a start point 
at radius < 30 cm from the beam line and an end point at radius > 37.5 cm. The total measured 
radial track length is required to be more than 10 cm and the specific energy loss to be more 
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than twice that of a MIP, i.e. log^o dE/dx > 0.3? These criteria are loose enough to ensure that 
the combined track reconstruction and particle identification efficiency is high (measured to be 
98.8 ± 0.2%). In the hard selection, tracks are required to pass through both CJCl and CJC2. 
The number of hits in CJCl is required to be at least 20 with a minimum of 40 hits in both 
chambers. The total measured radial track length must be greater than 35 cm. In addition, the 
selected tracks must have a minimum of 75% of the total number of possible hits and a specific 
energy loss of more than 2.5 times that of a MIP, i.e. log]^o(d-£'/da;) > 0.4. 

2.3 Particle Identification. 

The mass, M, associated with each charged particle track is deduced from the track momentum, 
p, and the most probable specific energy loss, dE/dxo, which is determined using a Bayesian 
log-likelihood method [13, 14]. The value of dE/dxo is chosen for each track such that the 
likelihood function 

N 

logL = ^\ogP{dE/dxi\dE/dxo) (2) 

i 

is maximised. Here, P{dE /dxi\dE / dxo) is the probability that the i^^ measured value of 
dE/dx results from a particle with most probable specific energy loss dE/dxo which is treated 
as a variable in equation El This probability is computed from a parameterisation of the Lan- 
dau distribution. The ratio p/M is then obtained in an iterative way from dE/dxo, assuming 
that the particle is singly charged, using a parameterisation of the Bethe-Bloch formula for the 
restricted energy loss [15] which includes corrections for apparatus effects. The value of M is 
then calculated using the measured track momentum. The maximum measurable mass using 
this technique is beyond the limit set by the centre of mass energy. 

Figure[l](upper plot) shows the specific energy loss for positively charged tracks, determined 
in this way, plotted against the track momenta. Clear bands can be seen corresponding to pions, 
kaons, protons, deuterons and tritons. The lower plot shows the spectrum of masses assigned 
to the tracks using the procedure described above. The smooth curves show parabolic fits to the 
log^o distributions (i.e. Gaussians) in each of the different mass peaks. Some deviations from 
Gaussian behaviour in the tails of the distributions are observed. The resolution is SM/M ~ 7% 
as determined from the widths of the Gaussian fits. The log-likelihood method adopted here is 
found to have better mass resolution and leads to a more Gaussian-like distribution than the 
method used previously by HI [16]. The particle type for each track is identified as the mass of 
the closest known particle. 

2.4 Background Determination 

The distributions of the z vertex coordinate and the distance of closest approach (DCA) of 
the tracks to the beam line in the transverse plane are used to distinguish the tracks produced 
in photoproduction from those produced by interactions of the beams with residual gas in the 
beam pipe, termed beam-gas interactions, or from the secondary interactions of photoproduced 

^Throughout the paper dE/dx is given as the ratio of the specific energy loss of the track to that of a MIP. 
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particles in the material of the beam pipe or the detector, termed material background. The 
beam-gas backgrounds are labelled pG for proton-gas and eG for electron-gas interactions. 

Figure 121 shows the DCA distributions. Tracks from interactions of the beam particles con- 
tribute to the peaks at zero, whereas the material background gives rise to the observed smooth 
background distributions. The material background for positive tracks results mainly from sec- 
ondary interactions with the material of the beam pipe or the other material before the CJC sen- 
sitive region. This can be seen in the upper plots of figure |2l where it leads to an approximately 
flat background with small peaks at DCA ~ 4 cm. For negative tracks, the material background 
is much smaller (see figure flower plots). It arises mainly from protons and deuterons which 
are back-scattered from the calorimeters into the CJC. These albedo particles, which lie outside 
the main peak in figure El are delayed due to their extra distance of travel by times of about 
6 ns relative to the arrival time of tracks coming directly from the photoproduction interaction 
vertex, as measured in the CJC. The selection of negative particles with masses greater than 
the proton mass is supplemented by the requirement that this delay be less than 4 ns, which is 
observed to reduce such backgrounds by about a factor of two. 

The number of particles corrected for the material background is obtained by subtracting 
the number in the sidebands of the DCA distribution (1.5 < |DCA| < 3.0 cm), normalised to 
the width of the selected region around the peak, from the total number of particles in the peak 
region. This width is chosen to be ±1.5 cm for p and p and ±0.5 cm for heavier particles for 
which the DCA resolution is better. 

The beam gas background is measured most accurately by dividing the data into four event 
samples depending on the presence or absence of energy in the forward part of the LAr or in the 
backward direction in the SpaCal. The event sample with both forward and backward energy 
(labelled Cn and comprising 89.4% of the total), arises dominantly from photoproduction (7^) 
which produces a roughly uniform distribution of energy in the apparatus for the VF^p range of 
this measurement. The event sample with forward but no backward energy (labelled Cio and 
comprising 8.5% of the total) arises dominantly from pG interactions, since protons interacting 
in the interaction region with a nearly stationary target produce mainly forward but little signif- 
icant backward energy. The event sample with backward but no forward energy (labelled Coi 
and comprising 2.0% of the total), is enriched in eG interactions for similar reasons. Only 0.1% 
of the events have neither forward nor backward calorimetric energy (sample Coo)- 

The data are separated into the ^p, pG and eG components in each sample, Cik, by studying 
the z vertex distributions (see figure |3l). The following procedure is adopted for the separation, 
the results of which are shown in table [l] for the hard selected tracks. First, the number of 
7P and pG (eG) tracks in the samples Ciq (Cqi), for each particle type, is measured by fitting 
the z vertex distributions to the sum of a Gaussian and a linear background. The integral of 
the Gaussian is taken to be the number of particles from 7p interactions, A'^^o (^of )' while the 
linear background determines the number of pG (eG) tracks, A^f^ (^oF)' '^he samples. Here 
the number of eG (pG) events in the Ciq (Coi) sample is neglected. The number of eG particles 
in the Cn sample is also negligible. When the distribution has no visible Gaussian shape (e.g. 
figure|3Jl), the value of N^^ is so low that a measurement is impossible (the dashes in tabled]). 
Second, the number of pG particles in each Cn sample, Nf^, is obtained by assuming that the 
probability that a true pG event has significant backward energy is small. Hence, the probability 
that a pG event appears in the Cn sample is approximately the probability that random noise 
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above threshold occurs in the SpaCal calorimeter in a pG event and this should be independent 
of particle type. This probability is measured from the ratio of the numbers of events in the Cn 
and Cio samples for pure pG events, i.e. samples which should have little contribution from 
photoproduction. Deuterons with z vertex \z\ > 20 cm and cos6' > are used for this, as is an 
independent sample of events with two identified protons and no identified anti-proton, which 
also shows no significant Gaussian shape in the z vertex distribution. The two independent 
measurements give ratios which agree within errors, confirming the assumption, with a mean 
of 0.151 ± 0.045. The quantity iVff is obtained by multiplying N^^ by this ratio. Finally, the 
number of 7p particles in the Cn sample, iV^f , is obtained by subtracting Nf*^ from the total 
for this sample, Nn. The separation into components for the Cqo sample is made using the 
probabilities (~ 0.02, estimated from the numbers in table [B that the photoproduced events 
have zero energy in the forward and backward calorimeters. 

2.5 Track Efficiencies 

The apparatus is fully sensitive in the range 0.2 < pt/M < 0.7 and \y\ < 0.4. In order to derive 
cross sections, corrections for track efficiencies (labelled e) must be applied. These are shown 
in detail for anti-deuterons in table|2lin the measurement intervals of px/M. The efficiencies for 
anti-protons tend to be somewhat larger than those for anti-deuterons, since the soft selection is 
used and the secondary interaction cross section is smaller. The efficiencies listed in table|2lare 
defined as follows. 

• ^dE/dx is the correction for migrations across the limit logdE/dx > 0.4. This is assessed 
by studying the migrations of anti-protons, selected with the looser criterion log dE /dx > 
0.3. 

• e,^ is a correction for a region of inefficiency in the CJC which developed during the data 
taking. 

• ecut is a correction for the loss of events outside the mass and DCA windows. 

• represents the corrections for the losses due to interactions in the material of the appa- 
ratus. This is extracted from the data using soft selected, identified tracks. The number 
of such tracks which are observed to interact in the material between CJCl and CJC2 is 
extrapolated to account for the material between the interaction point and the CJCl. The 
corrections agree with estimates from known cross sections. 

• ehit represents the track reconstruction efficiency which is determined by measuring the 
fraction of soft selected deuterons which enter the hard selection sample. 

• etrig is the trigger efficiency which is determined by two independent methods for events 
containing anti-protons. These are assumed to have the same trigger efficiency as anti- 
deuterons. One method involves Monte Carlo studies and the other method uses compar- 
isons of the number of events from the main trigger with those found by an independent 
monitor trigger. The two methods give consistent results. 

• ttag is the positron tagger acceptance. 
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• e^ch is the correction for the loss of events due to the requirement that there be five or 
more tracks in each event. This is deduced by applying the known KNO scaling distribu- 
tion [17], using a sample of events containing anti-protons which has an observed mean 
track multiplicity which matches that for events containing anti-deuterons. 

• e^o represents the correction for the losses of anti-deuterons outside the defined track 
timing interval. 

• ^PhSp is an estimate of the fraction of anti-deuterons lost outside the measurement region 
defined by the limits of |?7| < 1.0 and \y\ < 0.4. 

The final row of table|2lgives the measured value of the differential cross section 

da N 

= (3) 

d(pT/M) A{pt/M)CF ^ ' 

where is the number of events in each measurement interval, A{pt/M), corrected for all the 
efficiencies given in tableEl C = 5.53 ± 0.11 pb~^ is the integrated luminosity and F = 0.0136 
is the virtual photon flux per incident positron (for the details of the calculation see [18]). The 
first error quoted for the differential cross sections is statistical while the second is systematic, 
where the latter arises from the uncertainties in the efficiencies. All cross sections are quoted at 
the bin centres. 



3 Results 

3.1 Search for Heavy Particles 

The observed particle mass spectra are shown in figure IH The dashed curves show the mate- 
rial backgrounds deduced from the sideband subtraction method described in section l23l The 
mainly pG sample, Cio, contains 6 tracks which have reconstructed masses of more than 3 GeV. 
They each have specific energy losses which are approximately twice that expected for a MIR 
Visual inspection shows that these are overlapping tracks, probably due to relativistic particles 
which have been merged by the pattern recognition software, since the two tracks become visi- 
ble at the ends of their trajectories. No such tracks are seen in the dominantly photoproduction 
sample, Cn, which contains many more events. Since they are dominated by background, the 
Cio, Coi and Coo samples are omitted from the search for heavy particles in photoproduction. 

The observed deuterons and tritons (see tableland figure |H) are dominantly from the ma- 
terial background and upper limits on their photoproduction cross sections are derived from the 
observed numbers of events. These upper limits at the 95% confidence level, in the measure- 
ment range defined in section |231 are deduced to be 6.8 and 1.0 nb, respectively. No negative 
particles heavier than anti-deuterons and no positive particles heavier than tritons are observed. 
This allows an upper limit of 0.19 nb at the 95% confidence level to be set on the photoproduc- 
tion cross section for any such particle type in the same kinematic range. These cross section 
limits are derived assuming the same efficiencies as those for anti-deuterons given below. 
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3.2 The Anti-deuteron Cross Section 



A clear signal is seen in figure |4l consisting of a total of 45 anti-deuterons with an estimated 
material background of 1.0 ± 0.5. The inclusive cross sections are measured using the 35 
particles in the sensitive range defined in section 12.51 for all samples combined. Only two of 
the anti-deuterons are not in the Cn sample (both in Cio). This number is compatible with that 
expected from the probability that a photoproduction event has small backward energy (~ 2%). 

The total cross section for anti-deuteron production is found to be 2.7 ± 0.5 ± 0.2 nb in the 
kinematic range defined in section U31 by summing all the differential cross sections in the final 
row of table 121 The ratio of the number of anti-deuterons to anti-protons in the range \y\ < 0.4 
and 0.3 < pr/M < 0.7 is measured to be (5.0 ± 1.0 ± 0.5) ■ 10^1 The lower limit in pr/M 
is higher for the ratio than for the anti-deuteron cross section measurement in order to avoid the 
uncertainties associated with the large corrections to the anti-proton rates at low momentum. 

The inclusive anti-deuteron invariant cross section is given by: 

7 = M^E = (4) 

'd3(p/M) d3p 2n (pt/M) Ay d{pT/M) ^ 

where 7 = E/M, dcr/d(pT/M) are the differential cross sections given in the final row of 
table El and Ay = 0.8 is the rapidity range of the measurement. Figure |51 and table 13 show 
the measurements of this invariant cross section, normalised to the relevant total cross section, 
taken to be 165 ± 11/ib for photoproduction [18], as a function of px/M. Figure|51also shows 
the measured ratio of the corrected numbers of anti-deuterons to anti-protons versus px/M. 



3.3 Comparison with Other Measurements 

Measurements of the normalised invariant anti-deuteron cross section and the ratio of the anti- 
deuteron to anti-proton production rates performed by other high energy experiments in the 
central region are also shown in figure 13 The normalised cross section and ratio measurements 
obtained in pp collisions [5,6] and the photoproduction results described here are in good agree- 
ment, suggesting that the processes whereby anti-deuterons are formed are similar in pp and 7p 
interactions. The data on Au-Au collisions show a ratio which is slightly larger than that in pp 
collisions and photoproduction. However, the cross section for anti-deuteron production nor- 
malised to the total cross section in Au-Au collisions is over two orders of magnitude larger than 
that in photoproduction oxpp collisions, reflecting the more copious production of anti-nucleons 
in heavy ion collisions. 

The coalescence parameter B2 is derived by rearranging equation[l]to obtain 



Here, Mp and Mj are the masses of the p and d and a, ap and crj are the total interaction 
cross section and the partial cross sections for p and d production, respectively, as defined 
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in equation [T] with A = 2 for d production. The quantity, R, is the measured ratio of the 
number of anti-deuterons to anti-protons from direct production in each bin of pt/M, corrected 
for anti-protons formed remotely from the source by weak decays. In the coalescence model, 
only anti-nucleons produced directly from the source can form anti-deuterons. The number of 
p from direct production is taken to be 78 ± 8% of the number observed. This is estimated 
using the PYTHIA Monte Carlo [19], the accuracy being determined by the uncertainty in the 
strangeness suppression factor in this model. Hence, there is an overall theoretical uncertainty 
of about 20% in the determination of B2. The weak decay correction is somewhat larger for 
heavy ion collisions [7] which is thought to be due to enhanced strangeness production [20]. 

Figure Inland table |3l show the parameter B2, computed according to equation |5] from the 
data presented here. This quantity is calculated for the ISR data in [5,6], and is presented in 
figure |6l as a function of pt/M. The average value in photoproduction is determined from the 
data presented here to be B2 = 0.010 ± 0.002 ± 0.001 ± 0.002, where the first uncertainty is 
statistical, the second systematic from the errors in the efficiencies and the third the theoretical 
error from the weak decay correction. Figure|7]shows the mean value of B2 in photoproduction, 
compared with the mean B2 values obtained from the ISR pp data and the data of a variety 
of other experiments as a function of centre-of-mass energy. The value of B2 measured in 
photoproduction at W-yp = 200 GeV is similar in magnitude to the values deduced at lower 
centre-of-mass energies in pp and pA interactions (labelled "elementary" in figure|7l). However, 
this value of B2 is over an order of magnitude larger than that observed in Au-Au collisions 
at RHIC at a similar nucleon-nucleon centre-of-mass energy. Comparison of the heavy ion 
data with the data from more elementary targets shows that this discrepancy grows with centre- 
of-mass energy (see figure ITJ. To illustrate the difference between light and heavy colliding 
particles, the heavy ion data in figurelTJare restricted to very heavy ions. They are also restricted, 
at centre-of-mass energy below 50 GeV, to measurements of inclusive deuteron and proton 
production to avoid threshold effects in anti-deuteron production. The Bevelac data, at which 
energy the A dependence is weak, are the Ne-Au measurements of [21]. The AGS data are the 
Au-Pt measurements of the E886 experiment [22], the SPS data are the Pb-Pb measurements of 
NA44 [23] and NA52 [24] and the RHIC data are from the Au-Au measurements of the STAR 
Collaboration [7]. The "elementary" data are the pA data of [21, 22, 25], the pp data at the 
ISR [5, 6] and the photoproduction data presented here. 

In the coalescence model for heavy ion collisions, the parameter B2, varies inversely with 
both the volume of the fireball at thermal freeze-out and the absolute rate of anti-nucleon pro- 
duction when both the spatial and momentum dependence are included in the model [2] . The d 
to p ratio is observed in pp and 7J9 collisions to be close to that in Au-Au collisions, yet there is a 
much larger anti-nucleon production rate, and a smaller value of B2 in Au-Au collisions. These 
facts can be reconciled in the coalescence model if the size of the fireball at thermal freeze-out 
in pp and jp collisions is much smaller than that in Au-Au collisions. 



4 Conclusions 

A search for heavy charged particles is made in photoproduction at HERA and anti-deuterons 
are observed at (W^p) = 200 GeV. Upper limits at the 95% confidence level on the production 
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cross sections for any type of positive particles heavier than tritons or negative particles heavier 
than anti-deuterons are set at 0.19 nb in the kinematic range \y\ < 0.4 and 0.2 < px/M < 0.7. 
The total cross section for anti-deuteron photoproduction is measured to be 2.7 ± 0.5 ± 0.2 nb 
in the same kinematic range and the ratio of the number of anti-deuterons to anti-protons is 
measured to be (5.0 ± 1.0 ± 0.5) ■ lO"'' in the range \y\ < 0.4 and 0.3 < pr/M < 0.7. The 
transverse momentum dependence of the normalised invariant cross section for anti-deuteron 
production is found to be compatible with that measured in central pp interactions at a centre- 
of-mass energy of 53 GeV. The production rate per event of anti-deuterons in photoproduction 
is found to be over two orders of magnitude less than that observed in Au-Au collisions at 
RHIC, although the ratio of anti-deuterons to anti-protons is only slightly smaller. The coales- 
cence model parameter B2 is extracted in photoproduction and shown to be similar to that 
deduced from central high energy pp data and lower energy pA data. Averaging over the 
measurement range of px/M, B2 is measured in photoproduction at W^p = 200 GeV to be 
0.010 ± 0.002(stat.) ± 0.001 (sys.) ± 0.002(theory) GeVl This value is much larger than that 
observed in Au-Au collisions at RHIC at a similar nucleon-nucleon centre-of-mass energy. This 
difference between heavy ion and elementary particle collisions is reduced as the centre-of-mass 
energy decreases. These observations can be understood within the framework of the coales- 
cence model if the interaction volume at thermal freeze-out in 7p and pp collisions is much 
smaller than that in Au-Au collisions at a centre-of-mass energy of 200 GeV. 
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Figure 1: Upper plot - the observed specific ionisation energy loss, dE/dx, (normalised to that 
from a MEP), obtained from the log-likelihood method, versus track momentum (in GeV) for a 
sample of positively charged tracks from the hard selection (before application of the specific 

energy loss cut, see text). The smooth curves show the expected mean specific energy loss for 
the different particle species. Lower plot - the spectrum of masses (in GeV) for log^Q dE /dx > 
0.4, deduced as described in the text. The smooth curves represent Gaussian fits to each peak. 
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Figure 2: The distribution of the distance of closest approach (DC A) of the track to the event 
vertex in the plane transverse to the beam direction for p, d, p and d candidates with momenta 
larger than 0.5 GeV from the hard selection and before any track timing cuts. The peak at zero 
from beam induced events sits on the material background which is much smaller for p, d than 
for p, d. 
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Figure 3: The z vertex distributions for identified protons (a,b), deuterons (c,d), anti-protons 
(e,f) and anti-deuterons (g,h) in the hard selected samples after correction for the material back- 
grounds. The Cii samples (figures a,c,e,g) are dominated by photoproduction events and the 
Cio sample (figures b,d,f) are dominated by pG events. The solid curves (figures a,b,c,d) show 
the fits of a Gaussian distribution (expected from photoproduction) and a linear background (ex- 
pected from pG interactions). The lack of background and the relatively small number of events 
in (f) show that there is little p production in pG events. The number of d events in (h) is com- 
patible with that expected from the "fp contamination of Cio sample. The lack of background 
in (a) and (e) shows that there are few pG events in the Cn sample. The dashed histograms 
represent the background after subtraction of the fitted (b) or calculated (c,f) Gaussians. 
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Figure 4: The mass spectra (with M in GeV) for positive and negative particles in the hard 
selected sample for the mainly photoproduction (Cn) and mainly pG (Cio) samples. The dashed 
lines show the material background deduced from the sideband subtraction method. 
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Figure 5: a) The measured invariant cross section (normalised to the relevant total cross section) 
for inclusive d production, compared with the pp and Au-Au data, b) The measured d to p 
production ratio as a function of px/M (solid points) compared with pp data from the ISR at 
\/i = 53GeV in the central region [5, 6] and RHIC data on Au-Au collisions [7]. The inner 
error bars on the HI data indicate the statistical and the outer the total uncertainties. 
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Figure 6: The px/M dependence of the parameter B2 for 7p (this experiment), pp [5, 6] and 
Au-Au interactions [7]. The values of B2 for the pp data are deduced from the measured cross 
sections [5,6]. There is a theoretical uncertainty of about 20% (not shown) in the determination 
of B2 for the photoproduction and pp data, which arises from the calculation of the fraction 
of anti-protons produced directly (see text). The inner error bars on the HI data indicate the 
statistical and the outer the total uncertainties. 
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Figure 7: The dependence of B2 on centre-of-mass energy for very heavy ion collisions (open 
squares) and interactions of more elementary particles (closed circles). The heavy ion data are 
the Ne-Au data at the Bevelac [21], the E886 Au-Pt data at the AGS [22], the Pb-Pb data of 
NA44 [23] and NA52 [24] at the SPS and the Au-Au data of STAR [7] at RHIC. The "elemen- 
tary" data are the pA data of [21, 22, 25], the pp data at the ISR [5, 6] and the photoproduction 
data presented here. In the latter the error bar represents the total uncertainty i.e. the sum in 
quadrature of the systematic and statistical uncertainties. 
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Table 1 : The numbers (N^^'^ ) of photoproduction (7p) and beam gas (pG, eG) hard se- 
lected tracks (log^Q dE/dx > 0.4) in the four event samples (see text) for the different particle 
types. The right hand column gives the observed total number of identified tracks in each 
category after subtraction of the material background. The quoted errors are the statistical un- 
certainties. 
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Table 2: The observed number of events, track and event efficiencies and the differential cross 
section for J production (see section l231) . 
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Table 3: The measured values of the invariant d production cross sections, d to p ratios and 
coalescence parameter, B2. The first error is the statistical and the second error the systematic 
uncertainty. 
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